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Abstract 



We investigate the behavior of mixed 3 He- 4 He droplets on alkali surfaces at 
zero temperature, within the frame of Finite Range Density Functional theory. 
The properties of one single 3 He atom on 4 HeAf 4 droplets on different alkali 
surfaces are addressed, and the energetics and structure of 4 HeAr 4 + 3 He7v 3 
systems on Cs surfaces, for nanoscopic 4 He drops, are analyzed through the 
solutions of the mean field equations for varying number N% of 3 He atoms. 
We discuss the size effects on the single particle spectrum of 3 He atoms and 
on the shapes of both helium distributions. 

PACS 67.60.-g, 67.70.+n,61.46.+w 
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I. INTRODUCTION 



The physics of wetting by quantum fluids, particularly 4 He, 3 He and isotopic mixtures, 
received considerable attention from both experimentalists and theoreticians along the last 
decade, with experiments mostly pointing at measuring adsorption isotherms, determining 
the interfacial surface tensions, establishing the wetting temperature and constructing phase 
diagrams for the mixed systems. A geometrical, visible parameter that characterizes the 
wetting behavior of a liquid-substrate combination at temperature T is the contact angle of 
a macroscopic drop or wedge on a planar surface. The contact angle of liquid 4 He on Cs has 
been measured by at least three different groups, 1-4 and although the reported values do 
not coincide -a result which among other effects, could be attributed to the preparation of 
the surface- they are fully consistent with the helium inability to wet this weak adsorber. 5 
By contrast, 3 He has been theoretically shown to be a universal wetting agent 6 , a prediction 
which was verified shortly after, as isotherms of 3 He on Cs were measured that display neat 
prewetting jumps down to temperatures around 0.2 K. 7 

It is well-known that 3 He atoms admixed into bulk liquid 4 He or 4 He films at low T's 
populate a twodimensional (2D) homogeneous layer of Andreev states on the free surface 8 
(see also Ref. 9 and Refs. therein for a recent review). In this case, the Andreev state 
originates in the presence of a broad surface at the 4 He liquid-vapor interface and in the 
mutual interaction between the isotopes, and corresponds to the ground state (gs) of a dis- 
crete spectrum of states which become progressively localized towards the interior of the 
film. This structure appears as well in density profiles of mixed clusters. 10 ' 11 It was earlier 
proposed 12 ' 13 that if 4 He lies on a weak adsorber, Andreev-like states could be expected at 
the liquid-substrate interface. Theoretical anticipations of the wetting behavior of mixtures 
rely heavily on the expected reduction of the liquid-vapor surface tension of the mixture 
with respect to that of pure 4 He, due to the Fermi pressure of the Andreev surface layer. 14 
This effect, actually observed in Ref. 15, reduces the contact angle at a given temperature, 
thus permitting lower wetting temperatures as the 3 He concentration increases. On these 
grounds, Pettersen and Saam 14 predicted reentrant wetting of helium mixtures on Cs, a 
phenomenon measured shortly after by Ketola et a/. 16 The complete phase diagram for wet- 
ting of helium mixtures on alkali metal substrates derived in Ref. 17 offers a wide gallery 
of phenomena which includes prewetting, isotopic separation, triple point dewetting and 
A-transitions in the solution, in the full range of 3 He concentrations. Subsequent experimen- 
tal work discovered a rich pattern of wetting by mixtures, which includes, i.e., dewetting 
transitions near coexistence. 18 More recently, detailed measurements of the contact angle for 
dilute helium mixtures have been reported; 19 an analysis of the data showed that the large 
values of these angles are consistent with the presence of single-particle (sp) states of the 
3 He atoms, together with ripplons, at the liquid-solid interface. 

The first calculation of density profiles of finite droplets of 4 He on Cs, at zero temperature, 
was forwarded by Ancilotto et al. within the framework of Finite Range Density Functional 
Theory (FRDF) 20 , using the FRDF of Ref. 21. More recently, we have presented calculations 
of nanoscale mixed 3 He- 4 He droplets on Cs, 22 (hereafter, referred to as I) using a previously 
derived FRDF for mixtures. 10 ' 23 An important prediction in this paper is the existence of 
edge states, the lowest-lying 3 He sp bound states, which are essentially onedimensional (ID) 
and localize around the contact line of the 4 He drop -in contrast to the 2D surface Andreev 
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states. 

The purpose of this work is to examine in more detail the behavior of mixed helium 
droplets on alkali surfaces, at T=0. For this sake, in Sec. II we shortly review the current 
FRDF formulation, and in Sec. Ill we analyze the energetics and structure of a single 3 He 
atom on 4 HeAr 4 droplets on alkali planar surfaces, for varying number N± of 4 He atoms. In 
Sec. IV we discuss the solutions of the mean field equations for mixed 4 HeAr 4 + 3 HeAr 3 drops 
on Cs with varying number N 3 of 3 He atoms. Emphasis here lies on the features of the sp 
spectrum of 3 He atoms predicted by the solution of the Kohn-Sham (KS) equations derived 
from the FRDF, as well as on the shapes of both helium distributions. Our conclusions and 
perspectives are summarized in Sec. V. 



II. DENSITY FUNCTIONAL DESCRIPTION OF PURE AND MIXED HELIUM 
CLUSTERS ON ADSORBING SUBSTRATES 



The FRDF for helium mixtures adopted in this work is the same as in Ref. 23, of the 
general form 

E = j dv [U 3 (p 3 , r 3 , Pi) + W 4 (p 4 , r 4 ) + U u (p 3 , p 4 )] , (1) 

where Pi(r), Tj(r) are respectively the particle and the kinetic energy densities for i — 3, 4. 
The detailed form of the FRDF and values of the force coefficients have been given in Ref. 10, 
with the only changes reported in Ref. 23 which consist, on the one hand, in the neglect of 
the nonlocal gradient correction to the kinetic energy term in Ref. 21, and on the other hand, 
in the choice of the suppressed Lennard- Jones core as proposed earlier. 24 More specifically, 
for 3 He- 3 He, 3 He- 4 He and 4 He- 4 He interactions, the Lennard- Jones screened interaction is 
written as 



<Ju/r) 12 - {(Tu/rf if r > hi 



V L j{r) = < L J (2) 

{ V {r/hiP if r < hi , 

with e 44 = e 33 = e 34 = 10.22 K, <r 44 = <7 33 = 2.556 A, <r 34 = 2.580 A, hard-core radii 
K= 2.359665 A, h 3 = 2.356415 A and h 34 = 2.374775 A, and with V the value of the 
corresponding 6-12 potential at r = hi. These values have been fixed so that the volume 
integrals of the interactions Vj_,j coincide with the original ones in Refs. 10, 21. 
In the current geometry, for axially symmetric droplets, we have: 

p 3 (r,z) = 2 Y,\*m(r,z,ip)\ 2 (3) 

nl 

T 3 (r,z) = 2j2\Wni(r,z,tp)\ 2 (4) 

n I 

= , (5) 

4 p 4 (r,z) 



where 



e il(p 

*nj (r, z,(p) = ijira (r, z) — = (6) 
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denotes the sp wave functions (wf 's) of the 3 He atoms, and / denotes the projection of the 
sp orbital angular momentum on the z axis. Since we shall only address spin saturated 
systems, the sp energy levels are either twofold (I = 0) or fourfold (/ 7^ 0) degenerate. 

The Euler-Lagrange (EL) equations arising from functional differentiation of the density 
functional in Eq. (1) give rise to an integrodifferential coupled set 



h 2 

V 2 + V 4 (p 3 ,r 3 p 4 ) 



2m. 



y/pl = ha, sfpi (7) 



2m 3 



- W n * + V 3 (p 3 , r 3 , Pi) V n i = e nl (8) 



In Eq. (7), p 4 is the chemical potential that enforces particle number conservation in the 
4 He drop, and expression (8) represents the Kohn-Sham (KS) equations for 3 He atoms in 
the presence of their mutual interaction, their coupling to the 4 He cluster and the potential 
V s (z) generated by the alkali substrate filling the the z < half-space. For z > 0, the mean 
fields Vi{r, z) include V s (z) chosen as the Chizmeshya-Cole-Zaremba (CCZ) potentials. 25 The 
energies Si, i= 3, 4, of one helium atom in the substrate, obtained with the CCZ potentials, 
are shown in Table I. 

Eqs. (7) and (8) have been discretized using 7-point formulae and solved on a 2D (r, z) 
mesh. We have used sufficiently large boxes with spatial steps Ar = Az = h 4 /12 ~ 0.197A. 
As indicated in some detail in Ref. 26, we have employed an imaginary time method to 
find the solution of these equations written as imaginary time (r) diffusion equations. After 
every r-step, the 4 He density is normalized to N 4 , whereas the 3 He wf 's are orthonormalized 
using a Gram-Schmidt scheme. To start the iteration procedure, we have used the halved 
density of a 4 He27v 4 cluster obtained form a spherically symmetric FRDF code, and random 
numbers to build the 3 He wf's ip(r,z). This leaves little room to introduce any bias in the 
final results. 



III. ONE 3 HE ATOM ON 4 HE CLUSTERS SPREADED ON ALKAI SURFACES 

The first problem to consider is the sp spectrum for one single 3 He atom in the field of a 
4 HeAr 4 cluster on an alkali surface. As an example, we solve the corresponding Schrodinger 
equation for the 3 He atom for 4 He drops with N 4 = 20 and 100 on a Cs surface. In Fig. 1 
we show the sp level scheme e n \ in terms of squared angular momentum /, up to £3 = -3.13 
K (cf. Tab. I). For comparison, the lowest panel displays the sp spectrum obtained in I for 
a cluster with N 4 = 1000. Localized states with energy higher than -3.13 K are meaningless, 
as they would be an artifact of the calculation, carried out in a large but finite box. Indeed, 
on any substrate, 3 He atoms with energies higher than the corresponding £ 3 , would prefer 
to leave the 4 He neighborhood and occupy the lowest lying sp state of the alkali surface 
potential. One can see from Fig. 1 that whereas there are many 3 He sp states below S 3 for 
N 4 = 1000 (notice that only the lowest lying ones are shown), their number decreases with 
decreasing N 4 . In particular, for N 4 = 20 only the states with n — 1 and I — to 3 (s,p, d, 
and / states in spectroscopic notation), and the 2s state are bound to the 4 He droplet. 

In I we have shown that, for large N 4 values, states (nl) with n — 1 and / = 0, 1, 2, . . . 
are distributed into a rotational band. This means, on the one hand, that their sp energies 
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lie on a perfect straight line as functions of I 2 as seen in the bottom panel of Fig. 1, and 
on the other hand, that their probability distributions \ipu(r, z)\ 2 are sensibly identical to 
\ip w (r, z)\ 2 . Depending on N 4 , the same may happen, to some extent, to (2/) states with 
/ = 0, 1,2, . . ., (cf. Fig. 1). For N 4 = 20, the small size of the 4 He host cluster does not 
allow the energy levels to group into rotational bands, even for gs-based states ipu. This 
is due to the fact that the wf 's ip n o lie close enough to the z axis to experience centrifugal 
distortion at any nonvanishing /. The rotational character of the gs band en is recovered for 
iV 4 = 100, as depicted in the middle panel of Fig. 1. 

Contour plots of the probability densities [ip n0 (x, z)] 2 for n = 1 to 3, on the (x, z) plane, 
together with those of the density p 4 (x,z), are shown in Fig. 2 for N 4 = 100. It is clear 
that the gs -?/>io is localized on the circular contact line, revealing once again the edge state 
reported in I, where similar plots were presented for N 4 = 1000. In that work, we also 
showed that for n > 1, the probability densities display several fringes on the surface of the 
4 He cluster; a similar pattern occurs for N 4 = 100. This is not the case if N 4 is as low as 20, 
since although the edge state is still present, the small host immediately pulls out the 3 He 
probability density into unbound configurations (cf. upper panel in Fig. 1). We show in 
Fig. 3 contour plots of the probability densities [ipu(x, z)} 2 for I — 0, 1, and 2 (hereafter all 
the contour plots figures are drawn on the y = plane), together with those of the density 
p*{x,z), for N 4 = 20. 

Effective masses m* , defined 6 as the state averages of the parametrized local prefactor 
of the 3 He kinetic energy in the density functional Eq. (1) 

J-=[ dr ^i{r,z)] 2 

m* nl J m* 3 [p 4 (r,z),p 3 (r,z)} 

are displayed in Table II. The rotational character of the gs band en if N 4 = 100, allows 
one to fit these sp energies to a law 

£ "= £M+ 2Sk (10) 

with £10 = -4.52 K and regression unity. From m\ = 1.18 m 3 , we obtain a gs radius R\o = 
15.5 A, which sensibly coincides with the geometrical radius of the droplet at a 4 He density 
of about 10~ 2 A , as viewed in the bottom panel of Fig. 2. 

The appearance of the 3 He level structure built on these spreaded droplets is quite 
distinct. Within an strictly independent particle model, one may see from Fig. 1 that a 
fairly large amount of 3 He is needed before 21 states start being occupied. For iV 4 =100, 
only after filling the lh (/ = 5) state, i.e., N 3 = 22, the state 2s becomes occupied, and for 
A^ 4 =1000, only after filling the lk (I = 8) state, i.e., N 3 = 34. For this reason, the KS results 
discussed in Sec. IV refer to rather small 4 He droplets. Otherwise, the number of 3 He sp 
orbitals to compute, in order to see a sizeable effect due to the presence of this isotope, 
becomes prohibitive. 

It is also interesting to see how the number of ip n o localized bound states, with energies 
below the E 3 value in Table I, varies with N 4 and the alkali substrate. In particular, for Na 
and Li we have found only one single / = bound state, namely the edge state ipio, for N 4 
values from 20 to 3000. This is at variance with the results found for K, for which, similarly 
to Cs, the droplet with ^4=20 can only sustain one s state, whereas the drop with A^" 4 =100 
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and larger may sustain n — 1, 2 and 3 s states. We show in Fig 4 the contour plots of the 
probability densities [ip n o(x, z)] 2 for n — 1 to 3, together with those of the density p 4 (x,z), 
for N 4 = 3000 on K. As in the case of Cs 22 , the probability densities of the excited states 
display peaks on the upper surface of the drop. 

Contour plots of the probability density [tpio(x, z)] 2 , together with those of the density 
Pa(x, z), for N 4 = 3000 on Li and Na, are shown in Fig. 5. The fact that 4 He drops on Na 
and Li only support one bound s state, namely that surrounding the contact line, together 
with the rather weak / dependence of the sp energies in the ipu rotational band for large iV 4 , 
makes the 3 HeAr 3 rings hosted by these droplets, physical realizations of a neutral Luttinger 
liquid. 27 ' 28 

Figure 6 shows the energy E\ s of the edge state, as a function of N± for Cs, K, Na, and Li 
substrates, with ^(A^ = 0) = S3. It should be kept in mind that, as shown in Ref. 26, in 
the current FRDF description, Li, Na and K are wetted by 4 He, whereas Rb (not discussed 
here) and Cs are not. To establish a connection with the situation in films, for a given N 4 
droplet we define 

n 4 (r) = 2 7r J dzpi{r,z) (11) 

and in view of the analysis carried out in Ref. 26, it is safe to consider that 714(0) for the 
largest droplets A 4 = 3000 represents the prewetting coverage n c of a 4 He film on the given 
adsorber. This yields n c = 0.467 A" 2 for K, n c = 0.140 A' 2 for Na, and n c = 0.056 A -2 for 
Li. We have obtained the structure of films in the vicinity of n c , working out the first three 
3 He sp states as illustrated in Figs. 7 and 8. 

In Fig. 7 we plot the 4 He density as a function of the perpendicular distance z, for films 
on K, Na and Li, computed at the corresponding n c , together with the mean field for one 
3 He atom, the gs wave function tpo(z) and the first excited state wf tpi(z). In Fig. 8 we 
display the energies e n , n = to 2 (see also Fig. 7) for 4 He coverages around n c on the 
same three substrates, with the gs energy S3 of the 3 He atom on the substrate shown as a 
dotted line. The heavy dot indicates their corresponding n c . 29 As seen in Fig. 7, the wave 
functions follow the trend anticipated in Ref. 30: the gs of a 3 He impurity in a 4 He film 
remains mostly concentrated inside the Li monolayer, and is localized at the liquid-vapor 
interface if the film thickness is above one layer (K and Na). 

The energetics of the impurity shows an interesting evolution across the prewetting tran- 
sition. In the case of K, the weakest adsorber of this series, for the range of n 4 values 
here considered, the film spectrum is essentially constant with coverage, and for all three 
adsorbers the first excited state appears above the gs of the 3 He atom on the substrate. In 
contradistinction, at the prewetting density n c the energy of the Andreev state is below S3 
for K, and above it for Na and Li. We appreciate then that an important finite size effect 
for 4 He clusters on these later two adsorbers is the generation of the edge state for the 3 He 
impurity, which is favored with respect to binding to the substrate, thus permitting the 
existence of mixed drops of both helium isotopes on these adsorbers. It may also be inferred 
from Fig. 8 the different behaviour of adsorbers like K and weaker, as compared to Na 
adsorbers and stronger. Whereas in the former case, after filling the edge-like states, 3 He 
atoms start populating Andreev-like states covering the cup of the 4 He droplet, in the latter, 
once the edge-like states are filled, 3 He spreads on the alkali surface. 
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IV. MIXED CLUSTERS ON CS 



We now solve the coupled EL/KS equations (7) and (8) for N 4 = 20 and 100, and for 
several A3 values, on a Cs substrate. Fig. 9 displays the variation of the lowest KS sp 
energies e n i as functions of N 3 for the case A4 = 20, with heavy dots indicating the Fermi 
energy /chemical potential fi 3 of the given configuration. As seen in this figure, the chemical 
potential fi 3 is a -nonmonotonically- increasing function of the system size. This trend 
coincides with that encountered in 3 He, either pure 6 or dissolved in thick 4 He films. 12,13 
Figure 9 also shows that for N 4 =20, values of N 3 = 6, 10, 16, and 30 yield strong shell 
closures, and they are thus the magic numbers for the fermion component of a 4 He 20 + 3 HeAr 3 
droplet on Cs -these numbers depend on the alkali substrate and on A4-. Weaker shell 
closures also appear for A3 = 2, 20, 24, 34, and 38. 

Contour plots of the particle densities p 3 (x,z) and p 4 (x, z) are shown in Fig. 10 for 
A" 4 = 20, and in Fig. 11 for A" 4 = 100, and two values of A3. As we analyze these patterns, 
several observations may be put forward. On the one hand, we verify the spreading tendency 
of the density profiles p 3 with growing A3, as well as a tendency to cover the horizontal base 
and show some slight dilution inside the bulk of the host. It is also apparent from these 
figures that 3 He atoms in Andreev surface states do coat the 4 He cap. On the other hand, 
the shape of the 4 He droplet is rather insensitive to the presence of the intruder atoms, whose 
outwards spreading does not induce splashing of the 4 He cluster on the substrate surface as 
well. 

The behavior of the 3 He density should be interpreted together with examination of the 
KS level scheme in Fig. 9: an energetically based criterion for spreading of mixtures of 
(A3, A" 4 ) atoms on any substrate is given by the constraint p 3 (N 3c , N 4 ) = £ 3 , which defines a 
critical size N 3c (N 4 ) above which it is energetically more convenient, for an extra 3 He atom, 
to adhere to the flat substrate than to the 4 He drop. As seen in Fig. 9, for A" 4 =20, A3 C is 
about 44. 

Finally, a word of caution should be spoken with respect to the contact angle, for whose 
determination in large, macroscopic droplets, a procedure has been developed in Ref. 20. 
As discussed in I, due to the fact that nanoscopic 4 He droplet densities present large strat- 
ification near the substrate, the contact angle cannot be established by 'visual inspection' 
of the equidensity lines. Yet, even at a qualitative level of description, it is clear from the 
results here presented that addition of 3 He atoms affects the contact angle. To illustrate 
this assertion further, in Fig. 12 we show the contour plot of the total 4 He+ 3 He density 
for 4 He equidensity lines about 0.011 A , roughly half the bulk 4 He density, for A4 = 20 
and 100 and several N 3 values. Although the small size of these droplets prevents us from 
drawing quantitative conclusions, this figure indicates that the height of the cap changes 
less than the radius of its base, thus suggesting overall flattening of the combined density 
and a decrease of the contact angle, as expected according to both theory and experiment 
carried upon macroscopic samples. 

V. SUMMARY 

In our previous work here denoted as I, we have shown that in addition to the well known 
Andreev states, nanoscopic 4 He drops on Cs substrates host a new class of 3 He sp states, 



7 



edge states that generate a ID ring of 3 He atoms along the contact line. In this work we 
have carried a deeper investigation of size effects in mixed droplets on planar alkali surfacess, 
that enables us to assert that these edge states are a definite consequence of the particular 
combination of the confinement provoked by the adsorbing field, selfsaturation of the helium 
system and mutual interaction between isotopes. In fact, these localized states appear in all 
alkali substrates and for all 4 He clusters, although their ID nature becomes more evident 
either above nanoscopic host sizes or for the strongest confining potentials. We are able to 
compare the energy spectrum of single 3 He impurities in 4 He clusters with similar spectra on 
films, and found that while the evolution of the lowest lying sp states with film coverage near 
the prewetting transition may prevent solvation of a 3 He atom in a film on a strong adsorber 
such as Na or Li, its gs energy being unfavorable compared with binding to the substrate, 
the breaking of translational symmetry introduced by the finite size of clusters does change 
this trend, giving rise to an energetically favored edge state which remains localized around 
the drop. We also show that increasing N 3 /N 4 ratios enlarge the tendency of the mixed 
cluster to spread outwards on the surface; explicit KS calculations of the energy spectrum 
of added 3 He atoms and examination of the Fermi energy as a function of A3 support the 
quantummechanical interpretation for the onset of wetting by mixtures of helium fluids, as 
driven by the energetic benefit of abandoning the 4 He and adhering instead to the substrate. 
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TABLES 



TABLE I. Energies £3 and £4 (K) of one 3 He and one 4 He atom on different alkali substrates. 





Li 


Na 


K 


Rb 


Cs 


£3 


-9.81 


-6.41 


-3.74 


-3.31 


-3.13 


£4 


-10.70 


-7.08 


-4.20 


-3.72 


-3.53 



TABLE II. Effective masses m* (in units of m 3 ) for single 3 He atoms in He^ clusters on Cs. 



N 4 


20 


100 


500 


1000 


2000 


3000 


m* w 


1.13 


1.18 


1.20 


1.21 


1.22 


1.22 


m *20 


1.10 


1.26 


1.29 


1.30 


1.31 


1.31 


m *30 


1.05 


1.24 


1.30 


1.30 


1.30 


1.30 
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FIGURES 



FIG. 1. Single-particle energies e n i of one 3 He atom in a 4 HeAr 4 droplet on Cs, as functions of 
squared angular momentum projection, for N4 = 20, 100 and 1000. The dotted lines have been 
drawn to guide the eye. In each panel, from bottom to top, the states have principal quantum 
numbers n = 1,2,3. . ., and the dot-dashed line represents the gs energy £3 of one 3 He atom on 
the Cs substrate. 

FIG. 2. Contour plots of the probability densities [ip n o(x, z)] 2 for a single 3 He atom in a 4 Heioo 
droplet on Cs, together with the density P4(x, z) (bottom panel). In the case of 4 He, the equidensity 
lines correspond to 10~ 4 , 10" 3 , 2.5 x 10~ 3 , 5 x 10~ 3 , 7.5 x 10" 3 , 10" 2 , 2 x 10~ 2 , and 2.5 x 10~ 2 A~ 3 . 
In the case of 3 He, the equiprobability lines start at 10 _4 A 3 and increase in 10 _4 A 3 steps. 

FIG. 3. Contour plots of the probability densities [ipu(x, z)] 2 for a single 3 He atom in a 4 He2o 
droplet on Cs, together with the density p±(x, z) (bottom panel). The equidensity and equiproba- 
bility lines are as in Fig. 2. 

FIG. 4. Contour plots of the probability densities [ip n o(x, z)} 2 for a single 3 He atom in a 4 He3ooo 
droplet on K, together with the density p4(x,z) (bottom panel). The equidensity lines are as in 
Fig. 2, and the equiprobability lines start at 10~ 5 A and increase in 10" 5 A c steps. 

FIG. 5. Contour plots of the probability density [ipio(x, z)} 2 for a single 3 He atom in a 4 He3ooo 
droplet on Li, together with the density p^(x,z) (bottom panels), and on Na (top panels). In the 
case of 4 He, the equidensity lines are as in Fig. 2, and in the case of 3 He, the equiprobability lines 
correspond to 10" 8 , 10~ 7 , 10~ 6 , 10" 5 , 2 x 10~ 5 , 2.5 x 10~ 5 and 5 x 10" 5 A" 3 . 

FIG. 6. Energy £10 of the edge state (K) as a function of N4. The dotted lines have been drawn 
to guide the eye. The y axis is broken below ~ —7.5 K. 

FIG. 7. The 4 He density at the prewetting coverage n c on Na and Li substrates, and at 77,4 = 0.40 
A for K, together with the mean field V3 experienced by a 3 He impurity, and the gs wave function 
ipo(z) and first excited wf tpi(z) in that field. 

FIG. 8. The sp energies e n , n = 0, 1, and 2, for a 3 He impurity in 4 He films at coverages 
around n c on K, Na and Li. The dotted line indicates the gs energy £3 of one 3 He atom on the 
corresponding substrate, and the heavy dot points to the coverage at the prewetting jump. 29 

FIG. 9. The Kohn-Sham sp energies e n \ for 3 He7v 3 -|- 4 He2o mixed drops on Cs as functions of 
N3. The horizontal line indicates the energy £3 of one 3 He atom in the Cs adsorbing field. Heavy 
dots indicate the corresponding Fermi energy. The lines connecting the sp energies of (nl) states 
have been drawn to guide the eye. 
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FIG. 10. Contour plots of the density profiles p 3 (x, z) (top panels) and P4,(x, z) (bottom panels) 
of 4 He2o + 3 HeAT 3 mixed drops on Cs, for N 3 = 20 (left panels) and 40 (right panels). In the 
case of 4 He, the equidensity lines are as in Fig. 2, and in the case of 3 He, they correspond to 
10~ 4 ,2.5 x 10~ 4 ,5 x 10~ 4 ,7.5 x 10~ 4 , 10~ 3 , 2.5 x 10~ 3 ,5 x 10~ 3 ,7.5 x 10~ 3 , and 10~ 2 A~ 3 . 

FIG. 11. Same as Fig. 10 for N 4 = 100, and N 3 = 46 (left panels), and N 3 = 86 (right panels). 
The equidensity lines are as in Fig. 10. 

FIG. 12. Contour plots of the total 4 He+ 3 He density corresponding to 0.011 A 3 for the case 
iV 4 = 20, N 3 = 20 and 40 (top panel), and 7V 4 = 100, iV 3 = 46 and 86 (bottom panel). The dashed 
lines correspond to the pure 4 Hejv 4 drop. 
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